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We construct several new families of exactly and quasi-exactly solvable BCjv-type Calogero- 
Sutherland models with internal degrees of freedom. Our approach is based on the introduction 
of a new family of Dunkl operators of Bn type which, together with the original _Bjv-type Dunkl 
operators, are shown to preserve certain polynomial subspaces of finite dimension. We prove that a 
wide class of quadratic combinations involving these three sets of Dunkl operators always yields a 
spin Calogero-Sutherland model, which is (quasi-)exactly solvable by construction. We show that 
all the spin Calogero-Sutherland models obtainable within this framework can be expressed in a 
unified way in terms of a Weierstrass p function with suitable half-periods. This provides a natural 
spin counterpart of the well-known general formula for a scalar completely integrable potential 
[ of BCn type due to Olshanetsky and Perelomov. As an illustration of our method, we exactly 

. compute several energy levels and their corresponding wavefunctions of an elliptic quasi-exactly 

' solvable potential for two and three particles of spin 1/2. 
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I. INTRODUCTION 



I The completely integrable and exactly solvable models of Calogero and Sutherland describe a system of 
A'^ quantum particles in one dimension with long-range pairwise interaction. These models and their subsequent 
' generalizations (see |^ and references therein for a comprehensive review) have been extensively applied in many 
different fields of physical interest, such as fractional statistics and anyons ||, ^, quantum Hall liquids 0, Yang- 
1^ • Mills theories [|, |, and propagation of soliton waves . A significant effort has been devoted over the last decade 
, to the extension of scalar Calogero-Sutherland models to systems of particles with internal degrees of freedom or 
2 ■ "spin" In], ^ ^ ^ |8[ |o) . These models have attracted considerable interest due to their connection 
I ' with integrable spin chains of Haldane-Shastry type [^l], ^ through the "freezing trick" of Polychronakos |^ . 

The exactly solvable and integrable spin models introduced in |0, n4l generalize the original rational (Calogero) 
■ ' and trigonometric (Sutherland) scalar models, and are invariant with respect to the Weyl group of type Aj^. The 
exact solvability of both models can be established by relating the Hamiltonian to a quadratic combination of either 
the Dunkl |2j] or the Dunkl-Cherednik |^ operators oi An type, whose relevance in this context was first pointed 
out by Polychronakos We shall use the term "Dunkl operators" to collectively refer to this type of operators. 
^ ■ Up to the best of our knowledge, only two i?Ar-invariant spin Calogero-Sutherland models have been proposed so far 
\ in the literature, namely the rational and the trigonometric spin models constructed by Yamamoto in p6[ |. The exact 
solvability of the rational Yamamoto model was later proved in Ref. JlOf using the Dunkl operator formalism. The 
exact solvability of the trigonometric Yamamoto model will be proved in this paper. 

In a recent paper po|| the authors proposed a new systematic method for constructing spin Calogero-Sutherland 
models of type A^- One of the key ingredients of the method was the introduction of a new family of Dunkl- type 



operators which, together with the Dunkl operators defined in | |24| , 25 , preserve a certain polynomial module of 
finite dimension. It was shown that a wide class of quadratic combinations of all three types of Dunkl operators 
always yields a spin Calogero-Sutherland model. In this way all the previously known exactly solvable spin Calogero- 
Sutherland models oi A^ type are recovered and, what is more important, several new exactly and quasi-exactly 
solvable spin models are obtained. By quasi-exactly solvable (QES) we mean here that the Hamiltonian preserves 
a known finite-dimensional subspace of smooth functions, so that a finite subset of the spectrum can be computed 
algebraically; see j2^, |2^, ^ for further details. If the Hamiltonian leaves invariant an infinite increasing sequence of 
finite-dimensional subspaces, we shall say that the model is exactly solvable (ES). 



'Corresponding author. E-mail: artemio@eucmos.sim.ucm.es 

tOn leave of absence from Institute of Mathematics, 3 Tereshchenkivska St., 01601 Kyiv - 4 Ukraine 



2 



In this paper we extend the method of Ref. to construct new famihes of (Q)ES spin Calogero-Sutherland 
models of BC^ type. To this end, we define in Section II a new set of Dunkl operators of i?jv type leaving invariant 
a certain polynomial subspace of finite dimension, which is also preserved by the original Dunkl operators of Bm type 
introduced in |l9[| . In Section [II , we show that a suitable quadratic combination of all three types of Dunkl operators 
discussed in Section II can be mapped into a multi-parameter (Q)ES physical Hamiltonian with spin. This approach 
is a generalization of the construction used to prove the integrability of the An spin Calogero-Sutherland models, in 
which only a single set of Dunkl operators is involved. Our method is also related to the so-called hidden symmetry 
algebra approach to scalar A^-body QES models js^, where the Hamiltonian is expressed as a quadratic 

combination of the generators of a realization of s[{N + 1). We then show that the sets of Dunkl operators used in 
our construction are invariant under inversions and scale transformations. This property is exploited in Section IV to 
perform a complete classification of the i3Cjv-type (Q)ES spin Calogero-Sutherland models that can be constructed 
with the method described in this paper. The resulting potentials can be divided into nine inequivalent classes, out of 
which only two (the rational and trigonometric Yamamoto models) were previously known. In particular, we obtain 
four new families of elliptic QES spin Calogero-Sutherland models of BCn type. Section ^ is devoted to the discussion 
of the general structure of the potentials listed in Section We prove that all the potentials in the classification are 
expressible in a unified way in terms of a Weierstrass p function with suitable (sometimes infinite) half-periods. This 
provides a natural spin counterpart of Olshanetsky and Perelomov's formula for a general scalar potential related to 
the BCn root system. Finally, in Section VI we illustrate the method by exactly computing several energy levels and 
their corresponding eigenstates for an elliptic spin 1/2 potential in the two- and three-particle cases. 



II. Bjv-TYPE DUNKL OPERATORS 



In this section we introduce a new familiy of i?Ar-type Dunkl operators which will play a central role in our 
construction of new (Q)ES spin Calogero-Sutherland models. 

Let /(z) be an arbitrary function of z = (zi, . . . , zn) G K^. Consider the permutation operators Kij = Kji and 
the sign reversing operators i^i, whose action on the function / is given by 

iKijf){zi, . . . , Zi, . . . , Zj, . . . , zn) = f{zi, . . . , Zj, . . . , Zi, . . . , zn) , 
{Kif){zi,...,Zi,...,ZN) = /(zi,...,-Zi,...,zjv), 

where j, j = 1, . . . , TV. It follows that Kij and Ki verify the relations 

Kfj = 1, KijKjk = KikKi-j = KjkKik, Ki-jK^i ~ KkiKij, 

where the indices fc, I take distinct values in the range 1, . . . , A^. The operators Kij, Ki span the Weyl group of 
type Bn, also called the hyperoctahedral group. We shall also employ the customary notation Kij — KiKjKij. Let 
us consider the following set of Dunkl operators: 

Jt = ^^^-^^^ + AY. ^r^ (1 - ^^^) - E (1 - ^'^) 1 - ^'^41 - i-^TK,)^ (5) 
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where a, &, b' are nonzero real parameters, m is a nonnegative integer, and i = \, . . . , N . In Eqs. (|^)-(^, the symbol 
denotes summation in j with j = 1, . . . , i — 1, z -I- 1, . . . , iV. In general, any summation or product index without 
an explicit range will be understood in this paper to run from 1 to iV, unless otherwise constrained. It shall also be 
clear in each case whether a sum symbol with more than one index present denotes single or multiple summation. 

The operators in Eq. (^ have been used by Dunkl to construct a complete set of eigenvectors for Yamamoto's 
Bn rational spin model . The operators Ji were also introduced by Dunkl in Rcf. |jl^ . To the best of our knowledge, 
the operators jf' have not been considered previously in the literature. 
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The operators (|^)-(|^) obey the commutation relations 

2 

[Jf , Jf] = , [J°, J°] = ^ E + + Kjk - K,k - k,k) , (6) 



where e = ±,0, and the indices A: take distinct values in the range 1,...,7V. The operators (respectively 

i,j = 1, . . . , N, span a degenerate afhne 



Jj^), i — 1, . . . ,N, together with Kij and Ki, i,j = 1, . . . ,N, span a degenerate afhne Hecke algebra, see pa. The 



operators Jf do not commute, but since Jf + | J2j<ii-^ij + ^u) ~ § J2j>i(-^ij + ^u) can be shown that the 

latter operators, together with Kij and Ki, also define a degenerate affine Hecke algebra. 

It is well-known p9| that the operators J~ and Ji preserve the space Vn of polynomials in zi, . . . , zjv of degree at 
most n, for all n S N. Moreover, for any nonnegative integer n, the space TZn spanned by the monomials Yii ^1' with 
< k < n is also invariant under the action of both J~ and J°. Let us prove this assertion in the case of J?. Since 
{zi-^ — ^)TZn C TZn, it suffices to show that 



— (1 - A'y) n C 7^„ , and ^ (1 - ) [] 4" C Tin , (8) 



z, - z 



'3 J, ■^'l ^ ■^J 

for any pair of indices 1 < i j < N . For the first inclusion we note that 



'-l±^(l-K,^)]Jz[^^( ]lzi\z, + z,)iz,z,r'^(^^^^^^ ^ ' 

k \k^t,j ) 



\li Ij I \li Ij I 



ih - I,) ^ ^ 



Zj 



where 



l[ zi^Uz. + z,)iz.z,r'^<^'-'^hignik-l,) J2 
, fc#i,i / fc=o 



-1, P<0, 
sign(p) = <( 0, p = 0, 



(9) 



1, p>0. 

The resulting polynomial thus belongs to TZn- Indeed, it is a linear combination of monomials Y[k z]^ with 1^ = Ik for 
k ^ and li, Ij < max{li, Ij). Likewise, the second inclusion in ^ follows from the identity 



|:p^(i-^..)n4^= n zi^)iz.~-z,){z,z,r-^^''^^^s(i.,i,)^ — ' — 

\ k^i,j / fc=0 



(10) 



where 




We omit the analogous proof for the operators . Unlike the previous types of Dunkl operators, the operators 
m Eq. (|) do not preserve the polynomial spaces T'n and TZk with k ^ m. However, the space TZm is invariant under 
the action of . In fact, the operator 

zj-^-mz,-h'z,{\-{-\rK,), 
oz, 

preserves TZ„i, and, just as we did in the case of j", one can show that both 

(1 - K,, ) and (1 - K,, ) 



preserve TZn for any nonnegative integer n. 
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III. BCjv-TYPE SPIN MANY-BODY HAMILTONIANS 

In Section II we have shown that all three sets of B]\j-type Dunkl operators (||)-(^) preserve the finite-dimensional 
polynomial space TZm- In this section we shall use this fundamental property to construct several families of (Q)ES 
many-body Hamiltonians with internal degrees of freedom. 

Let 6 = span { |si, . . . , s^) | s., ^ -M, -M 1, . . . , M; M e 5N } be the Hilbert space of the particles' internal 
degrees of freedom or "spin" . We shall denote by Sij and Si, i, j — 1, . . . , N, the spin permutation and spin reversing 
operators, respectively, whose action on a spin state |si, . . . , sn) is defined by 

Sij \si, . . . , Si, . . . , Sj , . . . , S]y) — \si, . . . , Sj , . . . , Si, . . . , S]\f) , (11) 
Si\si, . . . ,Si,. . . , Sn) = |si, . . . , -Si, . . . , Sn) ■ 

The operators Sij and Si are represented in 6 by {2M + 1) ^-dimensional Hermitian matrices, and obey identities 
analogous to (|^). The notation Sij — SiSjSij shall also be used in what follows. 

We shall deal in this paper with a system of N identical fermions, so that the physical states are completely 
antisymmetric under permutations of the particles. A physical state ip must therefore satisfy AqiP = ip, where Aq 
is the antisymmetrisation operator defined by the relations Aq — Aq and n^jAo = — Aq, j > i ~ 1,...,N, with 
Ilij — KijSij. Since K^j = 1, the above relations are equivalent to Kij Aq — —Sij Aq, j > i — 1, . . . , N. For the lowest 
values of N, the antisymmetriser Aq is given by 

N = 2: Ao = i (1 - ni2) , 

N = 3: Ao = i (1 - - - + n^Uis + niaHaa) . 
6 

Our aim is to construct new (Q)ES Hamiltonians symmetric under the Weyl group of type Bn generated by 
the permutation operators 11^^ and the sign reversing operators KiSi. The corresponding algebraic eigenfunctions 
will be antisymmetric under a change of sign of both the spatial and spin variables of any particle, and therefore 
satisfy A^ — ip, where A is the projection on states antisymmetric under permutations and sign reversals. The total 
antisymmetriser A is determined by the relations A^ = A and 

K^JA = -S^JA, K^A^-S^A, J>l = l,...,N. (12) 

It may be easily shown that 

^^^(\{{1~KA)\Aq. 

^ i ' 

Following closely the procedure outlined in ||20[| , we shall consider a quadratic combination of the Dunkl opera- 
tors (|)-(|) of the form 

-H* = ^ (c++(J+)2 + coo( Jf)^ + c.-{J-? + CO Ji) , (13) 

i 

where c++,coo,c , cq are arbitrary real constants such that c^^ -I- Cqq + c?__ ^ 0. The second-order differential- 
difference operator (|l^) possesses the following remarkable properties. First, it is a quasi-exactly solvable operator, 
since it leaves invariant the polynomial space TZm- In particular, if c+_|_ = the operator H* preserves TZn (and Vn) 
for any nonnegative integer n, and is therefore exactly solvable. Secondly, H* commutes with Kij, Ki, Sij, and Si 
for alH, j = 1, . . . ,N. This follows immediately from the commutation relations (Q). Note that none of the terms 

Et-^'^]' Ei-^"^}' E^ 

i i i 

commute with Ki, and for that reason they have not been included in the definition of H* . We have also discarded 
the term YliiJt ^ ^i} because it differs from 2^ - [(J,f)^ + [h' — h)jf\ by a constant operator. 

Since H* preserves the polynomial module T^.^, commutes with A, and acts trivially on 6, the module 

= A(7^,„®6) (14) 
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is also invariant under H* . It follows from Eqs. (|T^) that the action of the operators Kij and Ki on the module TZm 
coincides with that of the spin operators —Sij and ~Si, respectively. Therefore, the differential operator H obtained 
from H* by the formal substitutions Kij —Sij, Ki ^ —Si, i,j — 1, . . . ,N, also preserves the module TZm- For the 



same reason, if the coefficient C++ in Eq. (13) vanishes, the operator H leaves the modules TZn and Vn = ^{'Pn (8) S) 
invariant for any non-negative integer n. Using the formulae (|^)-(^) in the Appendix for the squares of the Dunkl 
operators, we get the following explicit expression for the gauge spin Hamiltonian H: 

ZiP{Zi) 



i 

1 + s.,, 



z^j t ^ / 



(15) 



where 

P{z) = c++z^ + CooZ^ + c , 

Q{z) = 2c++(l - m - 6' + 2a(l - 7V))z3 + (cq + coo(l - m + 2a(l - A^)))z + 
i?(z) = c++TO(m - 1 + 26')z^ , 

—r + T2 I ^ ^4 - (^u- + S,j)^S,k + iifcjj + e 



26c__ 



(16) 



C = Coo 



?^ ( E' [4 - (^u- + + ~S^k)\ + 6 E(l - '^»^^) ) + ? E(2 + + ^y) 



Nmco 



Hereafter, the symbol Yl'ij k denotes summation in k with i j ^ k ^ i. 

One of the main ingredients of our method is the fact that the gauge spin Hamiltonian H can be reduced to a 
physical spin Hamiltonian 



(17) 



where T^(x) is a Hermitian matrix- valued function, by a suitable change of variables z = Ci'^)^ x — (xi, . . . ,xpf) and 
a gauge transformation with a scalar function /z(x), namely 



H-H\ 



(18) 



We emphasize that in general there is no (matrix or scalar) gauge factor and change of coordinates redu cing a g iven 
matrix second-order differential operator in N variables to a physical Hamiltonian of the form ([l7| ) ; see [ |33| , |34| and 
references therein for more details. The quadratic combination H* has precisely been chosen so that such a gauge 
factor and change of variables can be easily found for H. For instance, we have omitted the otherwise valid term 
^JJj^, J~] because it involves first-order derivatives with matrix- valued coefficients, which are usually very difficult 
to gauge away. The gauge factor jj, and change of variables z = C(x) in Eq. (^8|) are respectively given by 



= exp(5:/ 



my^) 



dyAl[{zf-z]rl[P{z,)-i. 



i<j 



and 



Xi= C ^{Zi) = 



dy 



The physical spin potential V reads 

a + Sij 



V^a 



E 



a + Si^ 



{Zi-Z.jf {Zi + Zj)"' 



i = l,. 



(19) 



(20) 



{zi + Zjf{l + S,j) + {zi - Zj)2(l -f Sij) 



b'c++ zf{l + (-1)™5,) + (1 + Sz) + W{z,) 



(21) 



C, 



6 



where 



P" 



4P 



P' 



3P' 



and 



C = aN{N -1) 



Co 



Coo 



(a(2iV-l) + 3(m-l)) 



-c, 



c = -c++(^2a'^{N - 1){2N - 1) + 4a(iV - !)(&' + m - 1) + m(26' + m - 1)) 



(22) 



(23) 



Note that the change of variables (po|), and hence the potential VYx), are defined up to an arbitrary translation in 
each coordinate Xi, i = 1, . . . , N. The hermiticity of the potential ( pi[ ) is a consequence of the Hermitian character of 
the spin operators Sij and 5'^. 

The invariance of the module TZm under the gauge spin Hamiltonian H and Eq. (18) imply that the finite-dimensional 
module 



Mr 



/i(x)A(7^„,(C(x)) ® 
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(24) 



is invariant under the physical spin Hamiltonian H. Therefore, any quadratic combination H* of the form (|l^) leads 
to a (quasi-)exactly solvable spin many-body potential (pTl)-(|23|) (provided of course that the module Aim is not 
trivial). In particular, if the coefficient c_|__|_ vanishes, the spin Hamiltonian H with potential ( ^l| ) is exactly solvable, 

since it leaves invariant the infinite chains of finite-dimensional modules Mn and Mn = ^J■{^) ^(Vn (C(x)) (8) &j , ti £ N . 

Our goal is to obtain a complete classification of the (Q)ES spin potentials of the form (pT| )-(p3[). The key observation 
used to perform this classification is the fact that different gauge spin Hamiltonians H may yield the same physical 
potential. This follows from the form invariance of the linear spaces span{ Jj^(z), J°(z), J^{z)}, i = 1, . . . , N , under 
projective (gauge) and scale transformations, given respectively by 



J^z) - J|(w) 



n 



,7V, 



±,0, 



(25) 



and 



Xz-i 



J/(z) ^ J|(w) = J^z), 



±,0, 



where A 7^ is real or purely imaginary. Indeed, we get 



jr(w) = -j+(w)| 



J«(w) = -Jl'(w), J+(w)=-J-(w)| 



b^b' 



(26) 



for the projective transformations (^5|), and J|(w) = A '^J/(w) for the scale transformations (|26|). This implies that 

the resulting quadratic combination H* is still of the form (p^, with (in general) different coefficients C++, cqo, c , 

and Co. Using these transformations, we can reduce the polynomial P{z) in ( |l5| ) to one of the following seven canonical 
forms: 



1, 

±vz^ , 
±iy{l + z'^), 

where 1/ > 0, < < 1, and < 9 < 7r/4. 



1. 
2. 
3. 
4. 



5. 
6. 
7. 



Ke'''-2')(e-''''-z'), 
±i^{l~z^){l~Pz^), 
z^){l-k^ + k^z^)^ 



(27) 



IV. CLASSIFICATION OF QES SPIN CALOGERO-SUTHERLAND MODELS 

We present in this section the complete classification of all the (Q)ES spin Calogero-Sutherland models that can 
be constructed applying the procedure described in the previous sections. To further simplify the classification, we 
note that the scaling (cee,co) ^-^ (Acee,Aco) induces the mapping 



(28) 
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of the corresponding potentials. For this reason, in Cases 2-7 we shall only list the potential for a suitably chosen 
value of the parameter v. Note furthermore that in Cases 2-4 and 6, once the potential has been computed for a 
positive value vq of the parameter v, its counterpart for the opposite value v = —vq can be immediately obtained 
using (^), namely 

F(x ; ~vq, Co) = -Viiyi ; vq, -cq) . 

For the models constructed to be symmetric under the Weyl group of type spanned by the operators KijSij 
and KiSi, 1 < i < j < N, the change of variables z = ^(x) should be an odd function of x, since only in this case 
X 1^ —X corresponds to z i-^ — z. In all cases except the second one, this has essentially the effect of fixing the 
arbitrary constants on which the change of variables (|20| ) depends. For example, in Case 7 with ly ~ A the change of 
variables is of the form Zi = ± C'n{2xi + | fc) = ± cn(2a;,; + ^i), 1 < i < N . Imposing that Zi be an odd function of 
Xi for all i and using the identity 

cn(2x,; + + cn(-2xj + ^j) = v v 



1 — sn^ sii^{2xi) ' 
we obtain the condition 

C^^{2h-1)K, keZ, l<i<N, 
where K = K{k) is the complete elliptic integral of the first kind 

Kik) 



V^l - fc2 sin^ ( 



Since cn{2xi~K+2liK) = (—1)'' cn{2xi—K), symmetry under exchange of the particles requires that k be independent 
of i, so that Zi — zt cn{2xi — K) for all i — 1, . . . , N. Taking into account that both the potential V and the gauge 
function fi are even functions of z by Eqs. (]l9|)-(p^), we see that the change of variables in this case can be taken as 
z, = cn(2xj - K), 1 < i < N. 

In the classification that follows, we have routinely discarded constant operators of the form 

= 70 + 71 51 ^» + ^2 ^ ^-^J + + S^j) + 74^'(5u + S,j){S.,k + S,k) , 7» e M. (29) 

i i<j i<j ij,k 

This is justified, since the operator Vq commutes with A (it actually commutes with Kij, Sij, Ki, and Si for \ < i < 
j < N) and therefore preserves the spaces A4„ and Nn for all n. 

All the potentials in the classification presented below are singular on the hyperplanes = Xj, 1 < « < j < A^, 
where they diverge as (xi — Xj)~'^. In some cases there may be other singular hyperplanes, near which the potential 
behaves as the inverse squared distance to the hyperplane. We shall accordingly choose as domain of the functions in 
the Hilbert space of the system a maximal open subset X of the open set 

xn < xn~i < ■ ■ ■ < X2 < Xi (30) 

containing no singularities of the potential. In all cases except Case 2b, we shall take as boundary conditions defining 
the eigenfunctions of H their square integrability on the region X and their vanishing on the boundary dX of X 
faster than the square root of the distance to the boundary. Since the algebraic eigenfunctions that we shall construct 
are in all cases regular inside X, when this set is bounded the square integrability of the algebraic eigenfunctions on 
X is an automatic consequence of their vanishing on dX. In Case 2b, the potential is regular and periodic in each 
coordinate in an unbounded domain. Therefore the square integrability of the eigenfunctions should be replaced by 
a Bloch-type boundary condition in this case. 

For each of the potentials in the classification, we shall list the domain chosen for its eigenfunctions and the 
restrictions imposed by the boundary conditions discussed above on the parameters on which the potential depends. 
In particular, the singularity of the potential at Xi — Xj, 1 < i < j < N , forces the parameter a to be greater than 1/2. 
Similarly, in all cases except for the second one the potential is also singular on the hyperplanes Xi = 0, 1 < i < N, 

-+(5 

and the vanishing of the algebraic eigenfunctions on these hyperplanes as |a;i| ^ with 5 > requires that b > 1/2. 
The conditions 

1 , 1 

^^^2' ^>2 
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shall therefore be understood to hold in all cases. For similar reasons, in Cases 4b, 5, and 6b we must also have 
The potential in each case will be expressed as 

i 

where the last term, which does not contain the spin operators Sij and Si, can be viewed as the contribution of a 
scalar external field. 

We shall use in the rest of this section the convenient abbreviations 



xf- =Xi±Xj, a = a{N - 1) + - {b + b' + m) . 



Case 1. P{z) = 1. 

Change of variables: z = x. 
Gauge factor. 



Scalar external potential: 



i<j i 



U{x)=co^x'^. (32) 

Spin potential: 

yspin(x) = 2a ^ [(x-.)-' (a + S^j) + (xt)-^ (a + Si^)] + b ^ xr^ {b + S,) . (33) 

i<j i 

Parameters: iv = —^cq > ■ 
Domain: < xn < ■ ■ ■ < xi. 

Case 2a. P(z) = 4:z'^. 
Change of variables: z = e^^. 

The most general change of variables in this case is Zi = Ae^^^* , 'i- < i < N, but the following formulas are independent 
of the choice of sign in the exponent and the value of the constant A. 

Gauge factor. 

Mx) = n[sinh(2a;r.)]". (34) 

i<j 

Scalar external potential: U (x) = . 
Spin potential: 

V'spiii(x) = 2a ^ l^sinh"^ a;~- (a + Sij) - cosh"^ x^j (a + Sij) . (35) 

i<j 

Parameters: cq = Am . 
Domain: xn < ■ ■ ■ < xi. 

Case 2b. P{z) = -4:z'^. 
Change of variables: z = e^'^. 

Again, the most general change of variables is Zi = Ae^^'^', 1 < i < -^j but the following formulas do not change 



9 



Mx) = n[sin(2a;r.)]V (36) 



when this is taken into account. 
Gauge factor. 



Scalar external potential: U (x) = . 
Spin potential: 

T4pin(x) = 2a ^ [sin-2 xr. (a + Sij ) + cos"^ a;^. (a + 5y ) . (37) 

Parameters: cq = —Am . 
Domain: xn < ■ ■ < x\ < xn + 

Both potentials in this case are invariant under a simultaneous translation of all the particles' coordinates. The choice 

Co = ±4m, which simplifies the form of the gauge factor, amounts to fixing the center of mass energy of the system. 
Note also that the potentials in this case do not possess symmetry, due to the fact that the change of variables 
cannot be made an odd function of x for any choice of the arbitrary constants. In fact, the sign change Zk ^ —zu 
corresponds to the translation Xk i— > a;/s +i7r/2 or Xk i— > a;fe +7r/2, which (as any overall translation) leaves the potential 
invariant. The potentials in this case are therefore best interpreted as ^Ar-type potentials depending both on spin 
permutation and sign reversing operators. 

For the hyperbolic potential 2a, none of the algebraic formal eigenfunctions are true eigenfunctions, since they are 
not square integrable on their domain. On the other hand, the algebraic eigenfunctions of the periodic potential 2b 
are clearly periodic in each coordinate and regular on their domain, and thus qualify as true eigenfunctions. 

Case 3a. P(z) = 4 (1 + z^). 
Change of variables: z = sinh(2a;). 
Gauge factor. 

m(x) = n [sinh(2x-.)sinh(2x+.)]" [J [sinh(2xi)]''[cosh(2a;,)]'^ . (38) 

i<j i 

Scalar external potential: 

U{x) = -Af3{(3 - 1) cosh-^(2a;) . (39) 

Spin potential: 

14pin(x) = 2a ^ |^( sinh"^ a;^. - cosh~^ xf-) (a + S^) + ( sinh"^ xf- - cosh~^ x^-) (a + Sij) 

(40) 

+ 46 ^ sinh-2(2a;i) (b+Si). 

i 

Parameters: /3 = ^ - (a{N -l) + b+^) <- {2a{N -l) + b + m) . 
8 V 2 / 

Domain: < xn < ■ ■ ■ < xi. 



Case 3b. P(z) -4 (1 + z^). 
Change of variables: z = i sin(2a;). 
Gauge factor. 

/x(x) = n [sin(2xr.)sin(24.)]" H [M^Xi)]''[cos{2xif . (41) 

i<j i 

Scalar external potential: 



U{x) =4/3(/3- l)cos-2(2a;). 



(42) 
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Spin potential: 

VspinCx) = 2aJ2[{ sin"^ + cos"^ (a + Sij) + ( sin-^ x± + cos"^ a;^.) (a + ^ij) + 4/; ^ sin-2(2a;i) (6 + 5i) . 

i<j i 

1 ^^^^ 

Parameters: p = - + a{N - 1) + 6 + ^) > - or /3 = 0. 

V 8 2/2 

Domain: < xn < ■ ■ • < xi < j, ii (3 > ^ and /3 ^ 1; < xat < • • • < < ^ — a;2, if /3 = 0, 1. 

Case 4a. P(z) = (1 - ^^^2 
Change of variables: z = tanh x. 
Gauge factor. 

m(x) = Yl (sinhxr. smhx±y JJe'^^^^^^a.,) (sinha;i)^(cosha;i)'''+'" . (44) 

i<j i 

Scalar external potential: 

U{x) = 20^ cosh(4a;) + 4/3(1 + 2a) cosh(2a;) . (45) 

Spin potential: 



Kpin(x) = 2a ^ j^sinh ^ x^^ (a + 5^ ) + sinh ^ a;+- (a + S'y)j 

+ 6 ^ sinh-^aji (6 + 5^) - 6' ^ cosh.-'^Xi {b' + (- 1)^5^) . (46) 

j i 

Parameters: /? = (cq + 2(6 - b')) < , or /3 = and a < 0. 
8 

Domain: < xjv < • • • < xi. 

Alternatively, we could have taken the change of variables as z = tanh {x — ^) = cothx. The gauge factor, external 
potential and spin potential become, respectively, 

m(x) = Yl ( sinhxr. sinhx+.)" JJ Q-0cosH2a^,) (cosh x^ )'' (sinh Xi)^'+™ , (47) 

i<j i 

U{x) = 2/?2 cosh(4x) - 4/3(1 + 2a) cosh(2x) , (48) 

and 

'54pin(x) = 2a ^ l^sinh"^ x,^- (a + Sij) + sinh~^ x^- (a + 5y)j 

- 6 ^ cosh-^Xi (6 + 5i) + 6' X] sinh-^Xi (b' + (-l)™^,) . (49) 



Case 4b. P{z) = -{1 - z'^f. 
Change of variables: z = i tanx. 
Gauge factor. 

/x(x) = Yl (sinx-- sinx+y JJe'^^^^^x,) (gin :ci )'' (cos x^ )*''+" • (50) 

i<j i 

Scalar external potential: 



U{x) = -20^ cos(4x) - 4/3(1 + 2a) cos(2x) . 



(51) 
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Spin potential: 

V"spin(x) = 2a ^ j^sin"^ x~j {a + 5^ ) + sin"^ x^l^ {a + 5^ )j 

+ b ^sin-^a;^ (b + Si) + b' ^cos-^a;^ {b' + (-l)'"5i) . (52) 



Parameters: = — „ (co + 2(6' — b)) . 
8 

Domain: < xn < ■ ■ ■ < xi < ^. 

The change of variable can also be taken as z = itan (-^ — f ) = i cot x. Since this is the result of applying an overall 
real translation to the particles' coordinates, we shall not list the corresponding formulas for the potential and gauge 
factor. 



2i0 2\/-2ie 2 



Cases. P{z) = {e^'^ -z^){e 
Change of variables: 



)■ 



sna; dna; 



cna; 



where the modulus of the elliptic functio ns is fc = cos 9. We shall also use in what follows the customary notation k' 
for the complementary modulus Vl — k'^- 

Gauge factor. 



/sna;..dna;..sna;+.dnx+-y j rfc 

^^"^ = n l-fc2,,2,-3,2,+ |/5-ctan [- cn(2x.)J } 



Scalar external potential: 



Spin potential: 



U{x)=4k'^dn-^{2x) 



X [sn(2a;i)]^[l + cn(2xi)] 



k0 



\{b'-b+m.) 



(3^ - a(a + 1) - -5- (1 + 2a) cn(2a;) 



[dn(2a;i)] (53) 



(54) 



14pin(x) = 2a ^ 



dn^ X 



1 + 



dn^ x% 



dn^ 



(a + Sij) + 



1 + 
sn^ xjj 



sr? Xi 
dn^ X. 



Parameters: ^ = ^ fco + 2(fc^ - k''^)ib - b')] . 

8kk \ I 

Domain: < xn < ■ ■ ■ < xi < K. 

An alternative form for the change of variables in this case is 

sn(a; — K) dn{x — K) cnx 
cn(a; — K) sn a; dn a; 

The resulting potential is obtained from the previous one by applying the overall real translation Xi ^ Xi — K, 
i = 1,...,N. Note also that, although Zj is singular at the zeros of cnx^, the algebraic eigcnfunctions satisfy the 
appropriate boundary condition on these hyperplanes on account of the inequality b' > 1/2 and the identity 

2 cn^ X 



(55) 



1 + cn(2a;) = 



1 — k'^ sn^ X 
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Case 6a. P{z) = 4{1 - z^){l - z^). 
Change of variables: z = sn(2a:). 

Here, as in the remaining cases, the Jacobian elhptic functions have modulus fc. 
Gauge factor. 

„ (sna;-.cna;--dna;-.sna;+.cna;^^dna;^^)" M^^ M/^' /rcA 

Mx) = ll ' \ ; _ ; 2, ^ ll[sn(2a;i)J [cn(2a;i)J [dn(2a;i)J • (56) 

Scalar external potential: 

U{x) = 4fc'2 [/?(/? - 1) cn-2(2a;) - - 1) dn-2(2a;)] . (57) 

Spin potential: 



V;pin(x) = 2a ^ 



cn^ x,^, dn^ x,-; ,. sn^ a;,t \ / cn^ x^! dn^ x,^- ,. sn' x. 



sn2xy cn2a;tdn xt.y \^ sn2a;^ cn^ a;^^. dn a;^ y 

+ 46 ^ sn-2(2xi) (& + S^) + Ak%' ^ sn2(2a;i) (&' + (-l)'"S'i) . (58) 

i i 

Parameters: 

/3 = -^(co+4(l + fc2)(6-6')) -a> i or /? = , = _L (cq + 4(1 + _ 5')) _ a . 

Domain: < < ■ ■ ■ < Xi < if /? > ^ and /? 1; < x^r < • • • < a;i < -fT — 0:2, if /3 = 0, 1. 

Case 6b. P{z) = -4(1 - _ ^2)_ 
Change of variables: 

. sn(2x) 
z = 1- 



cn{2x) 
Gauge factor. 

„ (snxy cnx,--dnxy snx+-cnx+-dnx+-)" Mfe'+"»rj/o M/^' /rm 

Mx)=ll ' ^ ' , _ „ .Aa ^ 11 [sn(2a;i)] [cn(2xi)] [dn(2x0] • (59) 

i<j [l-k^sn^x^jsn^xlj) ^ 

Scalar external potential: 

U{x) = A[l3{(3 - 1) k^ sn2(2x) - /?'(/?' - 1) k'^ dn-2(2x)] . (60) 

Spin potential: 



'i4pin(x) = 2a ^ 



+ . 2 + (« + -^ij) + 0+0+ + , 2 - (« + ^ij) 



dn^xt ) ysn^x+.cn^xf. dn^ x: 



+ 46 ^ sn-2(2a;i) (6 + 5^) + Ak'%' ^ cn-2(2a;i) (6' + (-1)'"^^) . (61) 



Parameters: 



l^=^{^o+ 4(1 + k'^){b' - b)) -a, (3' = -^{co+ 4(1 + k'^){b' -b))-a. 



Domain: < xn < ■ ■ ■ < xi < ^ 
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In spite of the singularity of Zi at the zeros of cn(2a;i), the vanishing of the gauge factor on these hyperplanes clearly 
implies that the algebraic eigenfunctions fulfill the appropriate boundary condition. 



Case 7. P(z) = 4(1 - z^)(fc'^ + fc^ z^) 
Change of variables: 



Gauge factor. 



z = cn(2a; - K) = k 



I sn(2a::) 
dn(2a;) 



m(x) = n 

Scalar external potential: 



Spin potential: 

X4pin(x) = 2a ^ 



sn cn x,^ dn a;,, sn xti cn xt^ dn x. 



— dn ^ 



(l — fc2 so? xfj] 



Yl [sn(2x,)]' [cn(2a;0]'' [dn(2a;,; 



1 b'+m 



U{x) = 4 [/?'(/?' - 1) s,u^{2x) + - 1) cn~2(2x)] . 



2 + 1 2 

sn X -' dn 



sn2 x,^. dn^ x^ 



(a + 5'y) + 



2 + 

cn^ a; 



^ sn2 0:+- dn^ 



■ 45 ^ sn-2(2a;,) (6 + 5^) - ^ dn"2(2a;,) (6' + (-1)"S'.) . 



sn^x^^ dn x^^- 



(62) 



(63) 



(a + ^y) 



(64) 



Parameters: 



j3^^['^ + he -k'^){b' ~h) + o\>\ or /3 = 0, = ^ + i(fc2_A:'2)(6'-6)-a. 



Domain: 0<xjv<---<xi<-y, if/3>i and /^t^I; Q < xn < ■ ■ ■ < xi < K — X2, if/3 = 0, 1. 



V. DISCUSSION 



The -BCAf-type potentials constructed in the previous section can be expressed in a unified way that we shall now 
describe. In the first place, apart from irrelevant constant operators of the form (p9|), the spin potential can be written 

as 

F,pin(x) = 2a ^ [(t;(x-.) + v{x+ + Pi)) (a + %) + («(x+) + v{x~^ + P^)){a + 5,,) 



i<3 



b H^'O + «(^. + Pi)) (b + S^) + b'J2 (^(^* + ^2) + vix, +Pi+ P2)) (&' + (-1)™5,), 



(65) 



where w is a (possibly degenerate) elliptic function, and Pi and P2 are suitably chosen primitive half-periods of v (see 
Table H). In particular, when one of the periods Pi of v goes to infinity, expressions like v{x + Pi) with x G K finite are 
defined as zero. In Case 1, both periods are infinite and w(x-|-Pi-|-P2) is also defined as zero. Furthermore, the constant 
K' = K'{k) in the elliptic Cases 5-7 is the complete elliptic integral of the first kind defined by K'{k) = K{k'). Using 
this notation, it is easy to verify that in the non-degenerate elliptic Cases 5-7 the scalar external potential U{x) can 
be written as 



U{x) = A(A - 1) v{x + iPi) -I- v{x - iPi) + A'(A' - 1) v{x + iPi + P2) + v{x ~ \Pi + P2) , (66) 

where X = f3 and A' ~ (3' in Cases 6-7, while A = —a + i/3 and A' = —a — i/3 in Case 5. Formula (^6|) holds also in 
Case 3, with A = /3 and A' = 0. In Cases 1 and 4 Eq. (|6^ ) cannot be directly applied, since in these cases all the terms 
in (pq) are either indeterminate or zero. However, the potentials in Cases 4a and 4b can be obtained from that of 
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Case 


v(x) 


Pi 


Pa 


1. 




oo 


oo 


3a. 


sinh^^ X 


m 

T 


oo 


3b. 


smT^ X 


TT 

2 


oo 


4a. 


sinh~^ X 


00 


m 
T 


4b. 


siri^^ X 


oo 


TV 

2 


5. 


sn^ X 


K + \K' 


K 


6a. 


cv? X dv? X 
sn^ X 


K 


iK' 
2 


6b. 


dr? X 
sn^ X cv? X 


iK' 


K 
2 


7. 


CV? X 

sn^ X dv? X 


K 





TABLE L Function v{x) and its primitive half-periods Pi (see Eq. (pq)) for each of the PCjv-type potentials in Section 



Case 5 in the limits — > and — > 7r/2, respectively. Likewise, applying the rescaling Xi t-^ vxi (i = I, . . . , N , v > 0) 
to the potential of type 3a or 3b one obtains the potential of type 1 by taking /3 = — a>/(4i/^) and letting v ^ 0. 

The function v{x) that determines the potential V^(x) in the elliptic Cases 5-7 according to Eqs. ( |65| ) and ( |66| ) can 
be expressed in a systematic way in terms of the Weierstrass function p{x; uji, 073) with primitive half-periods cui = K 
and CJ3 = iK' . Indeed, dropping inessential constant operators we have 

v{x) = e[p{x; wi, W3) + p{x + 2P2; wi, 1^3)] , (67) 

where P2 is the primitive half-period of v listed in Table |, and e — 1 for Cases 6-7 while e — 1/2 for Case 5 (the 
only case in which 2P2 = 2K is a period of p). Since in Cases 6-7 Pi and 2P2 are primitive half-periods of p, the 
well-known second-order modular transformation of the Weierstrass function Q applied to Eq. ( |67| ) leads to the 
equality 

v{x) = pix;Pi,P2), (68) 

where the primitive half-periods Pi and P2 are listed in Table ^ and we have dropped an irrelevant additive constant. 
Substituting Eq. (|6^) into Eqs. ( |65| ) and ( |66| ) and applying once again a modular transformation to the one-particle 
terms we readily obtain the following remarkable expression for the potential V^(x) in Cases 5-7: 

y(x) = 2a ^ [(p(x-.; Pi, P2) + p(4 + A; Pi, P2)) (a + 5,,) + (p(a;+-; Pi, P2) + p(Xy. + Pi; Pi, P2)) (a + ) 

+ 46 ^ p(2x,; Pi, 2P2)(6 5.) -f 46' ^ p(2x, -I- 2P2; Pi, 2P2)(6' + (-1)"^,) (69) 



4 ^ [a(A - 1) p(2x, + Pi; Pi, 2P2) + A' (A' - 1) p(2a;, + Pi + 2P2; Pi, 2P2) 



One of the main results in this paper is thus the fact that the potential (^9|) is QES provided that the ordered 
pair (Pi,P2) is chosen from Cases 5-7 in Table |. In fact, the remaining PC^v-type (Q)ES spin potentials listed 



in Section IV can be obtained from the potentials in Eqs. ( |69| ) by sending one or both of the half-periods of the 
Weierstrass function to infinity. This is of course reminiscent of the analogous property of the integrable scalar 
Calogero-Sutherland models associated to root systems ||] . 
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The potentials in Cases 1, 2, and 3 are ES for all values of the parameters. (In Case 3, the dependence of the 
parameter /3 on m through a can be absorbed in the coefficient Cq.) The potentials of type 4 are also ES for (3 — 0. 
The elliptic potentials in Cases 5-7 are always QES. 



All the potentials presented in Section IV are new, except for Cases 1 and 4. Case 1 is the rational BAr-type model 
introduced by Yamamoto [|l^ and studied by Dunkl jl^]. Case 4b for /3 = is Yamamoto's BAr-type trigonometric 
potential with Ai = —b (in the notation of Ref. jl^), and either X[ = —6' < —1/2 for m even or A'^ = 6' > 1/2 for m 
odd. Our results thus establish the exact solvability of the trigonometric Yamamoto model when \X[\ > 1/2. 

It should be noted that the method developed in this paper admits a number of straightforward generalizations. In 
the first place, the algebraic states could be chosen symmetric under sign reversals. The resulting Hamiltonians would 



coincide with the ones presented in Section [V with 5*^ replaced by —Si. In particular, if 6 = &' = one can obtain 
algebraic eigenfunctions of both types (symmetric and antisymmetric under sign reversals) for the same Hamiltonian. 
The construction can also be applied to a system of N identical bosons, just by replacing the antisymmetriser Aq by 
the projector on states symmetric under permutations of the particles. Choosing a system of fermions is motivated 
by the fact that the internal degrees of freedom can be naturally interpreted as the physical spin of the particles when 
M = 1/2. 



The procedure described in Section |I| relies on the algebraic identities analogous to (|^) satisfied by the spin 
operators Sij and Si, and not on the particular realization (p^). For instance, replacing the operators Sij by new 

operators Sij spanning one of the anyon-like realizations introduced by Basu-Mallick [T^ would yield further families 
of (Q)ES spin Calogero-Sutherland models. 

VI. EXACT SOLUTIONS FOR AN ELLIPTIC QES MODEL 

As an illustration of the procedure described in the previous sections, we shall now compute the algebraic sector of 
the spectrum for the elliptic QES potential of type 6a in Eqs. (^-(|8|) in the case of two and three particles of spin 
1/2 (N = 2, 3 and M = 1/2), for m = 1, 2, 3. Note that in the spin 1/2 case, the spin permutation and sign reversing 
operators Sij and Si can be expressed in terms of the usual one-particle SU(2) spin operators cr^ = {al,af, af) in the 
more familiar way 

Sij = 2(Ti ■ (Tj + Si — 2(j] . 



The operator H* corresponding to the potential (|57|)-(p8D reads 



N 



H*^-J2 (4fc^J+)^ - 4(1 + k')iJ^f + A{jrf + CO jf) + C* + E, , (70) 

i=l 

where C* is the constant operator obtained by replacing Sij by —Kij and Si by —Ki in the expression ( [T^ ) for C, 
and the scalar constant Eq is given by 

Ea =coN{a{N -l) + b' + m + ]^ - 2N{1 + k'^){2a[N - l)[2b' + m) + 2b' [b' + ni + I) + ni +'^{N - l)[2N - l)a^^ . 

Let us first consider the two-particle case [N — 2), for which the spin space S is spanned by the four spin states 
|±±) = |±i ±i). For m = 1 the polynomial module TZi is the one-dimensional space span{(/3i}, with 

^1 = (.1 - ^2)(|+ + > - I—)) + {Zl + Z2){\-+) | + -» . (71) 

Therefore, the spin state 

V'i(x) = ^/^l [sn2;i2Cna;+ dna;^2(|-h+) - | )) + sna;+ cnxj^a dnxj^j (|-+) - (72) 

1 — fc^ sn^ sn '- ^ ' J 

is an eigenfunction of the Hamiltonian of type 6a, where the gauge factor /i(x) is given in Eq. (^6|). The corresponding 
eigenvalue is Ei = Eq — Cq. 

If m = 2, the antisymmetrised polynomial module TZ2 is the three-dimensional space spanj^si, ip2, fs}, where (pi is 
given by Eq. ( [7l|) and 

^2 = {zf - Zl){\++) + I—) - I + - |- + » , 

ip3 = ZiZ2{zi - Z2){\++) - I )) + ZiZ2{zi+Z2){\ + -) - |-+)) . 
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The matrix of the gauge spin Hamiltonian H (or H*) in the basis </?2, fz] is given by 



/£;o - Co - 4(1 + P) 8(26+1) 

So - 2co 

\ 8(2&' + l)fc2 So - 3co - 4(1 + fc2)y 



whose eigenvalues are 



Si^3 = So - 2co - 4(1 + fc2) ^ A , S2 - So - 2co , 
where A = [co + 64fc^(26 + 1)(26' + 1)] • The corresponding physical wavefunctions are 



V'l,3(x) = 



1 — A;^ sn^ a;]^2 sn^ 
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sna;^2Cna;^2dna;j*2(coT A + 8A;^(26' + 1) sn(2xi) sn(2a;2) ) (|++) - | )) 



snXj^ cnX]^2 dna;j^2 



■(/'2(x) = A*(x) 



(co T A - 8fc2(26' + 1) sn(2xi) sn(2a;2)) (|-+) - |+-)) 
sn xTn cn xTn dn xTn sn xfn cn a;^'r, dn cc^, 



(74) 



(l — fc^ sn^ 



^<2) 



(l++> 



I + - 



-)) 



For m = 3 the antisymmetrised polynomial module TZj, is spanned by the spin functions (/Si, . . . , cpg, where 

= {zl - z|)(l++> - 1—)) + {4 + ^l)(l-+> - 1+-)) , 

(^5 = ^1^2(2? - 2|)(l + + > + I — > + l + -> + 1- + )) , 

= zlzlizi - Z2){\++) - I )) + zlzl{zi + Z2)(|-+) - 1 + -)) . 

The matrix representing H — Eq in the basis {ipi, ■ ■ ■ , ipe} is 

/-co-8(l + fc2) 


8(2&' + 3)/c2 







8(26+1) 
-2co-8(l + fc2) 

-3co - 16(1 + /c^) 



V 



-8(26' + 1)P 








8(26' + l)fc2 



-8(4a + 26 + 3) 


-16a(l + /c2) 
~3co + 16a(l + /c2) 


-8(4a + 26' + 3)fc2 



-4co - 8(1 + A:2) 








8(26 + 3) 
-8(26+1) 


5co-8(l + fc2)y 



(75) 



Clearly, So — 2co — 8(1 + fc^) and Sq — 4co — 8(1 + k^) belong to the spectrum of H and hence of H, with corresponding 
eigenf unctions given respectively by V'2(x) in Eq. ([74|) and 



, , , , , sna;,, cnxi, dnxno sna;^, cna;^, dnx^, 

^5(X) = /i(x) _ 'I ^2 

(1 — sn^ sn^ X12) 



i(2xi) sn(2a;2) (|++) + |— ) + | + -) + |-+)) 



(76) 



The remaining algebraic levels are the roots of a fourth degree polynomial, whose expression is too long to display 
here. For instance, if a = 6 = 6' = 1, A:^ = 1/2, and cq ~ —14 (so that f3 = 0), the algebraic levels are approximately 
Si = -327.4, S2 = -288, S3 = -281.4, S4 = -262.2, S5 = -260, and Sg = -201.0. 

In the three-particle case, the spin space 6 is spanned by the eight states |±±±). If m ~ 1, the antisymmetrised 
space TZi is trivial. For m — 2, the antisymmetrised space TZ2 is spanned by the single state 

"Pi = Z12 Z23 (!+ + + > + I )) - Z^2 A-i ^23 (!+ + -> + I +» 

- 4l 43 ^23 (1 + — ) + !-++» + 42 ^13 43 (! + - + > + 1- + -)) ' (77) 

where zf^ — Zi ± Zj. Consequently, tpii'^) — /^(x)!/?! (z), with Zi — sn(2a;i), is an eigenfunction of H with eigenvalue 
So-3co-4(l + fc2). 
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When m = 3, a basis for TZ^ is given by the function (pi in Eq. (|7^) and 

(P2 = Zfs ^23 (^1 + ^2 + Z3) (!+++> - I >) + 2+ (zi + Z2 - Z3) (| h) - |+ + ->) 

+ ^12 ^213 i-^l +Z2 + Z3) (1 + ) - |- + + )) + Z+ Zfs 2^3 (^1 - Z2 + Z3) (| + -+) - | - + -)), 

V3 = Z^2 ^13 ^23 (^1^2 + Z1Z3 + Z2Z3) (I + ++) + I >) + ^13 ^23 (^^1^2 + ^1^3 + Z2Z3) (| h> + |+H )) 

+ ^12 ^1*3 ^23 (^1^2 + 2:1^3 - ^2^3) (|H ) + I h + )) - zf^ ^13 ^23 (^1^2 - Z1Z3 + Z2Z3) (|H h) + | \ )) 

ip4: = Z1Z2Z3 {z^2 ^13 ^23 (! + + + > - I )) + ^12 43 ^23 (!++-> " I + >) 

+ 42 43 ^23 (h + + > - I + — » + 42 H3 43 (|- + -> " ! + - + >)] • 

The matrix of H in the basis {</?!,..., <y94} reads 

/i;o - 3co - 16(1 + 8(4a + 26 + 3) \ 

i;o-4co + 8(l + P)(2a-l) 8(26+1) 

8(26' + l)fc2 £;o-5co + 8(l + fc2)(2a-l) 

8(4a + 26' + 3)fc2 - 6co - 16(1 + fc^)/ 

The eigenvalues of this matrix are 

Ex^3 ^Eo- 4co + 4(1 + k^){2a - 3) T A_ , E2,a ^ Eq ~ 5co + 4(1 + k^)(2a - 3) T A+ , 

where 



(78) 



V 



A+ = (co + 4(1 + fc2)(2a +l)y + 64fc2(26 + l)(4a + 26' + 3) 



1/2 



1/2 



A_ = (co - 4(1 + fc2)(2a + 1)) + 64fc2(26' + l)(4a + 26 + 3) 
The corresponding eigenfunctions are 

^i^3(x) = m(x) [(co - 4(1 + fc2)(2a + 1) T A_) (^i(z) + 8^^(26' + 1) ^3(z) 

V'2,4(x) = /x(x) (co + 4(1 + fc2)(2a + 1) + A+) ^2(2) + 8fc2(4a + 26' + 3) ^4(2) 

with Zi — sn(2a;i). 



Appendix 



In this Appendix we present a hst of identities satisfied by the Dunkl operators (||)-(H) used to compute the gauge 
spin Hamiltonian (111). 



z^ - z^ 



a 



-,2 ' 



Y.^j'f = E (^'. + (1 - + 2a(l - A^))z.9.,) + 4a 5] 9.. + aj] 



z2 - z2 



(Z, - Zj) 



ZiZj iVm^ 



{k,, - 1) + + ^ [4 - (X,, + x,,)(if,, + k^k)] + 6 E(i - ^^^iO 



(z, + Zj)2 ^ ' ■ 4 ' 12 



Y.^4f = E (^'. - + ^ - 1 + 2a(^ - l))z!d.^) + 4a 5] -2^ + a^ 



z^z^ 



Z^Z^ 
Zi z^ 



(z, + Zj)- 



(X., - 1) + 6' ^ ((-l)"i^. _ 1) + ^(^ _ 1 + 26') ^ zi , 



Nm 



The symbol ^- ^ means summation in z, j, A: with i ^ j ^ k ^ 
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